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Abstract. We derive the critical behavior of a CA traffic flow model using an 
order parameter breaking the symmetry of the jam-free phase. Random braking 
appears to be the symmetry-breaking field conjugate to the order parameter. For 
Vraax = 2, we determine the values of the critical exponents /3, 7 and 5 using an 
order-3 cluster approximation and computer simulations. These critical exponents 
satisfy a scaling relation, which can be derived assuming that the order parameter 
is a generalized homogeneous function of \p — Pc\ and p in the vicinity of the phase 
transition point. 
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1 Introduction 



Cellular automaton (CA) models of traffic flow have attracted much interest 
since the publication of the Nagel-Schreckenberg (NS) model fl]]. The NS 
model is a probabilistic CA model of traffic flow on a one-lane highway. The 
road is represented by a lattice of L sites with periodic boundary conditions. 
Each site is either empty (in state e) or occupied by a car with velocity 
V = 0,1,2,... ,fmax (in state v). If di is the distance between car i and 
car i + 1 (cars are moving to the right), velocities are updated in parallel 
according to the following two subrules: 

Vi{t + 1/2) = mm{viit) + 1, di{t) - l,v^,,) (1) 
Vi{t + 1) = max(t>j(t + 1/2) — 1,0) with probability p, (2) 

where Vi{t) is the velocity of car i at time t; then, cars are moving according 
to the subrule: 

Xi{t+1) =x^it)+v,{t + l), (3) 

where Xi{t) is the position of car i at time t. The model contains three 
parameters: the maximum speed fmax, which is the same for all cars, the 
braking probability p, and the car density p. For a clear presentation of the 
various approximate techniques to calculate relevant physical quantities of 
the NS model compared with results obtained from computer simulations 
see §. 

If p = 0, the NS model is deterministic and the average velocity over the 
whole lattice is exactly given by 

{v) = min ^Wmax, ~ ^ ^) • (4) 

This expression shows that, below a critical car density pc = l/(fmax + 1), all 
cars move with a velocity equal to fmax, while above Pc, the average velocity 
is less than fmax- This transition from a free- moving regime to a congested 
regime is usually viewed as a second-order phase transition. 

If p 7^ 0, except in the case v^ax = 1, no exact expression for the average 
velocity has been obtained, and there is no consensus concerning the existence 
of the phase transition (see, for instance and 0] and references therein). 
One point that is particularly unclear is the definition of the order parameter. 
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If CA models of traffic flow exhibit second-order phase transitions, we 
have to understand the nature of the order parameter, show how it is related 
to symmetry-breaking, determine the symmetry-breaking field conjugate to 
the order parameter, define the analogue of the susceptibility, study the crit- 
ical behavior, and, possibly, find scaling laws. 

The purpose of this short note is to fulfill this program. In our discussion 
we shall consider a slightly simpler model of traffic flow. We shall assume 
that the acceleration, which is equal to 1 in the NS model, has the largest 
possible value (less or equal to Wmax) as in the Fukui-Ishibashi (FI) model 0, 
except that these authors apply random delays only to cars whose velocity 
is equal to fmax- That is, in our model, we just replace (0) by 

Vi{t + l/2) = mm{di{t) - 1, v^^y,). (5) 

2 Symmetry considerations 

Deterministic CA rules modeling traffic flow on one-lane highways are number- 
conserving (p = constant). Limit sets of number-conserving CA have, in most 
cases, a very simple structure and, these limit sets are reached after a number 
of time steps proportional to the lattice size |P, 0, H]. The limit set of our 
model is identical to the limit set of the NS or FI models. 

If p < pc, any configuration in the limit set consists of "perfect tiles" of 
■^^max + 1 cells as shown below 
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e 




e 


^max 



in a sea of cells in state e. In particular, it can be shown that the prob- 
ability to find in the limit set a sequence of f max + 1 empty sites is exactly 
given by 

P(ee-_e) = 1 - (t;max + l)p- 

This relation is a simple consequence of the following argument: If, in an 
empty road, you add Lp cars, in the limit set, for p < pc, we have 

LP{^e_^-e) = L - (w^ax + 

where L is the lattice size. 

If p > Pc, a configuration belonging to the limit set only consists of a 
mixture of tiles containing v + 1 cells of the type 
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where v = 0, 1, ■ ■ 


5 ^max- 
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= 0, 


1,2,- 



J is the velocities 



distribution, we have 




Note that (§) is a simple consequence of these relations. 

If we introduce random braking, then, even at low density, some tiles 
become defective, which causes the average velocity to be less than fmax- 
The random-braking parameter p can, therefore, be viewed as a symmetry- 
breaking field, and the order parameter, conjugate to that field is 



(6) 



This point of view implies that the phase transition characterized by m will 
be smeared out in the presence of random braking as in ferromagnetic systems 
placed in a magnetic field. 

From (Bl) and (H), it follows that, for p = 0, 



m 





P- Pc 

PPc 



if P < Pc, 
otherwise. 



(7) 



Hence, the critical exponent (3 is equal to 1. 



3 Approximate techniques 

To determine the other critical exponents, we have used local structure ap- 
proximations as described in [10] and numerical simulations. 

To construct a local structure approximation, it is more convenient to 
represent configurations of cars as binary sequences, where zeros represent 
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empty spaces and ones represent cars. Since for Vmax = 2 the speed of a car 
is determined by the states of, at most, two sites in front of it, the minimal 
block size to obtain nontrivial results is 3 (the site occupied by a car plus 
two sites in front of it). In what follows, we limit our attention to order-3 
local structure approximation. 

Using 3-block probabilities, we can write a set of equations describing the 
time evolution of these probabilities 

Pt+iihhh) = ^ w{b2b3h\aoaia2a3a4a5afy)Pt{aoaia2a3a4a5aQ), (8) 

^^£{0,1} 
i=0,l,... ,6 

where Pt{h2h3h4) is the probability of block &2&3&4 at time t, 
and W7(62&3&4|aoaia203a4a5a6) is the conditional probability that the rule 
maps the seven-block into the three-block 62^3^4- Letters 

a denote the states of lattice sites at time t, while h denote states at time 
t + 1, so that, for example, 03 is the state of site i = 3 at time t, and 63 is 
the state of the same site at time t + 1. Conditional probabilities w can be 
easily computed from the definition of the rule, although, since their number 
is quite large, we used a computer program for their determination. 

Equation (H) is exact. The approximation consists of expressing the seven- 
block probabilities in terms of three-block probabilities using Bayesian ex- 
tension. That is, 

Pt{aQaia2a3aia^aQ) = (9) 

Pt(aoaia2)Pt(QiQ2Q3)-Pt(Q2Q3a4)-Pt(Q3Q4a5)-Pf(Q4Q5a6) 
Pt{aia2)Pt{a2a3)Pt{a3a4)Pt{a4a5) 

where Pt{aiai+i) = Pt(ajaj+iO) + Pf(ajaj+il) for i = 1,... ,4. Equations 
d^) and @ define a dynamical system whose fixed point approximates three- 
block probabilities of the limit set of the CA rule. Due to the nonlinear nature 
of these equations it is not possible to find the fixed point analytically. We 
have determined it numerically. 

From the knowledge of three-block probabilities Pt{ ), the car den- 
sity is given by 

p = Pt{l) = Pi(lOO) + Pi(lOl) + Pi(llO) + Pi(lll), (10) 
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and the average velocity by 

^ 2(1 - p)Pt{m + pPtim + (1 - p)Pt{m ^^^^ 

This last result gives the expression of the order parameter m = Vmax ~ (v)- 
Note that p = is a constant of motion for equation (H) whose fixed 

point depends, therefore, upon the initial car density -Po(l) = -Po(lOO) + 

Po(lOl) + Po(llO) + Po(lll)- Consequently, using order-3 local structure 

approximation allows to determine the order parameter m as a function of 

car density, and the critical exponents. 

For most simulations, we used a lattice size equal to 1000 and our results 

are average of 1000 runs of 1000 iterations. For p = 0.0005 we took a lattice 

of 10000 sites and averaged 500 runs of 10000 iterations. 



4 Results 

The susceptibility Xp ^'^ constant p, defined by 

dm 

Xp = lim— , 12 

p^o op 

cannot be calculated exactly. So we determined its value using local structure 
approximation as described in the preceding section and computer simula- 
tions for t>max = 2. Figures l^a and |pD show the results obtained from order-3 
local structure approximation^ and simulations, respectively. In the limit 
p — > 0, the susceptibility diverges as (pc — p)~'^ for p < pc, and as (p — Pc)~'^ 
for p > pc- Using local structure approximation (Figures ^ and ^), we 
found 

7 = 0.91 ± 0.03 and 7' = 0.98 ± 0.03 
while simulations (Figures 0b and ||b) yield 

7 = 0.86 ±0.05 and 7' = 0.94 ± 0.05. 

Another exponent of interest is 6. It characterizes the behavior of m as 
a power of p for p = pc- Here again we have determined the value of 

^■^ "^(Pc,0) -m{pc,p) 
p^O p 

^If Vmax > 2, one should use (wmax + l)*''-order local structure approximation 



6 



using order-3 local structure approximation and simulations. Our results are 
represented in Figure ^, and the values of 6, obtained using order-3 local 
structure approximation and computer simulations, are, respectively, given 
by 

1/5 = 0.51 ±0.01 and 1/5 = 0.53 ± 0.02 

In is interesting to note that the values 

P = l, 7^1, 6^2 

obtained for the critical exponents are found in equilibrium statistical physics 
in the case of second-order phase transitions characterized by nonnegative 
order parameters above the upper critical dimensionality. 

Close to the phase transition point, critical exponents obey scaling rela- 
tions. If we assume that, in the vicinity of the critical point (p = PcP = 0), 
the order parameter m is a generalized homogeneous function of p — pc and 
p of the form 

where the function / is such that /(O) 7^ 0, then, differentiating / with 
respect to p and taking the limit p —* 0, we readily obtain 

7 = 7 = (5 - (14) 

which is verified by our numerical values obtained either using order-3 local 
structure approximation or computer simulations. Figures ^ and clearly 
confirm the existence of a universal scaling function. 

5 Conclusion 

In this short note, we have derived the critical behavior of a CA traffic 
flow model using an order parameter breaking the symmetry of the jam- 
free phase. Random braking, which is thought to be an essential ingredient 
of any CA traffic flow model, appears to be the symmetry-breaking field 
conjugate to the order parameter. For Vmax = 2, we have determined the 
critical exponents f5, 7 and 6 using order-3 local structure approximation 
and computer simulations. These critical exponents satisfy a scaling relation 
which can be derived assuming that the order parameter is a generalized 
homogeneous function of \p — Pc\ and p in the vicinity of the phase transition 
point. 
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Figure 1: Divergence of the susceptibility at the critical point using local 
structure approximation (a) and data obtained from simulations (b). Ex- 
cellent agreement of simulations and local structure approximation is clearly 
visible. 
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Local structure results 
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Figure 2: Plot of [m{p) — m{0)]/p as a function of |p — pd for P < Pc using 
local structure approximation (a) and data obtained from simulations (b). 
Straight line represents best fit for p = 0.0005 and 0.01 < pc — p < 0.1. 
We obtained 7 by performing such a fit for several different values of p and 
extrapolation to p = 0. 
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Local structure results 
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b) 

Figure 3: Plot of [m{p) — m{0)]/p as a function of |p — pd for P > Pc using 
local structure approximation (a) and data obtained from simulations (b). 
Straight line represents best fit for p = 0.0005 and 0.01 < p — pc < 0.1. 
We obtained 7' by performing such a fit for several different values of p and 
extrapolation to p = 0. 
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Local structure results 
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Figure 4: Order parameter m at the critical point p = pc as a function of p 
using local structure approximation (a) and data obtained from simulations 
(b). Slope of the best fit line at p = pc yields 1/S. 
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Figure 5: Illustration of the existence of universal scaling function for p < Pc 
using local structure approximation (a) and data obtained from simulations 
(b). 
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Figure 6: Illustration of the existence of universal scaling function for p > Pc 
using local structure approximation (a) and data obtained from simulations 
(b). 
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